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ABSTRACT. Results of tests to demonstrate the very high differential-photometric stability of CCD light 
sensors are presented. The measurements reported here demonstrate that in a controlled laboratory 
environment, a front-illuminated CCD can provide differential-photometric measurements with 
reproducible precision approaching one part in 10 2 * * 5 . Practical limitations to the precision of 
differential-photometric measurements with CCDs and implications for spacebome applications are 
discussed. 


1. INTRODUCTION 

Even though CCDs have now superseded other light sen- 
sors for most astronomical photometry, there is widespread 
uncertainty as to their attainable precision. Limitations pro- 
duced by external factors are sometimes thought to be intrin- 
sic limitations of the devices themselves, giving unnecessar- 
ily pessimistic predictions of what is possible with CCD 
sensors. Earlier work bearing on this subject includes papers 
by Buffington et al. (1990, 1991), Gilliland et al. (1988, 
1991), and Young et al. (1991). 

The measurements reported here demonstrate that under 
ideal circumstances, a front-illuminated CCD can provide 
differential-photometric measurements with reproducibility 
approaching one part in 10 5 . Although photometric precision 
for astronomy is normally much poorer due to factors such 
as the limited number of photons available and atmospheric 
scintillation, it is important to understand how to separate the 
limitations produced by the detector from limitations due to 
external and possibly correctable factors. 

The work reported here was prompted by a proposal to 
detect planets orbiting other stars by observing transits of the 
planets across the stellar disks. The transits would be de- 
tected by making precise measurements of photometric 
variations for a large number of stars over an extended pe- 
riod of time from a space platform, where many of the prob- 
lems seen by ground-based observers can be avoided. This 
paper describes the experimental tests of one CCD, then dis- 
cusses how the results were analyzed and how they may be 
interpreted. In addition, possible sources of photometric error 
not measured by this experiment are mentioned. 

2. EXPERIMENTAL DETAILS 

In order to demonstrate that precise measurements are 

possible with a CCD, images in the laboratory were taken of 

an artificial star-field over a ten-day period. The goal of the 
experiment was to show that one can detect small temporal 
changes in the light level of any one star when compared to 
the average of the other stars in the field. 

The experiment was carried out at the Lick Observatory 


CCD laboratory using a front-illuminated Reticon 512X512 
CCD with 27 micron pixels. The device was operated in 
multipinned-phase (MPP) mode at a stabilized temperature 
of about -110°C, and had a full well of approximately 
3X10 5 electrons per pixel. The output signal from the CCD 
was attenuated to produce about 40 electrons per DN, and 
signal levels of about 6000 DN were used to avoid saturating 
the CCD. A 16 bit analog-to-digital converter (ADC) was 
used for the measurements. The CCD had good charge- 
transfer characteristics for both large and small signals, al- 
though there were a few low-level charge traps scattered 
over the detector’s sensitive area. In order to reduce the ef- 
fects of photon shot-noise, 100 successive exposures were 
added together for each recorded image, resulting in approxi- 
mately 40X6000X 100— 24X 10 6 electrons per pixel in the 
final images. 

The apparatus consisted of the CCD mounted in a liquid- 
nitrogen-cooled dewar, with a simple illumination system to 
produce artificial star images. The CCD position was refer- 
enced to the front face of the dewar, rather than to the liquid- 
nitrogen cold surface, to reduce mechanical motion as the 
coolant evaporated. The CCD was illuminated with light 
from a single red LED that passed through two separated 
translucent plastic diffusers, then through a mask consisting 
of a metal plate drilled with holes of five different sizes. The 
mask was imaged onto the CCD by a Canon // 1.8 lens. The 
intensity of illumination across the field was uniform to 
within about 30%. In order to avoid jarring and other errors 
that might be caused by a mechanical shutter, the exposures 
were controlled by switching the LED on and off with a 
relay. The assembly was mounted with the optical axis ver- 
tical and placed in a large enclosure whose temperature was 
controlled to within about 0.1 °C. Metal rods were used to 
clamp the optical assembly to the dewar while cylinders of 
PVC with 1 in. thick walls were used as spacers for the 
optical setup. 

The liquid-nitrogen-cooled dewar had a hold time of only 
about 20 h, so regular refilling was required during the 10- 
day run reported here. As will be noted from the data pre- 
sented later, the LN refill was a major cause of mechanical 
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Fig. 1— ’ This shows the size and distribution of artificial stars imaged on the 
CCD. Only data from the marked medium-size “stars” were analyzed for 
this report. 

instability. In order to limit this problem, the fill funnel was 
carefully fastened in place and isolated as well as possible 
from the rest of the equipment. The temperature-dependent 
resistance of a carbon resistor mounted inside the dewar al- 
lowed a rough estimate of the LN level, to allow refilling 
without overflow. 

3. EXPERIMENTAL RESULTS 

Although a number of preliminary test runs were carried 
out in which various problems with the experimental setup 
were discovered and corrected, only data from the final 10- 
day run are presented here. Frames of 451 X460 pixels were 
taken repeatedly and the average of each succeeding 100 
frames was recorded as a single image; each recorded image 
represents about 35-min worth ot data. Data were recorded 
continuously over the 10-day period, producing 354 recorded 
images. In order to reach an arbitrary shot-noise limit of 1 
part in 10'^ it was necessary to add the data from about 400 
pixels, so “stars” of about that size were chosen for the 
analysis. The “star field” is shown in Fig. 1 and the cross 
section of a typical “star” used in the analysis is shown in 
Fig. 2. 

To correct for small unavoidable variations of the light 
intensity and amplifier gain, data from each star were divided 
by the average signal level from all the stars. Raw data and 
normalized data from a typical star are shown in Figs. 3 and 

4. Note the sudden changes when the dewar was refilled. The 

large spike in the normalized brightness following image 50 
(Fig. 4) corresponds to a single instance when the dewar was 
overfilled, resulting in spillage that produced thermally in- 
duced image motion. The motion of the centroid of the same 
star is shown in Fig. 5(a). Although the LED was driven by 
a voltage-regulated power supply, a 0.5% variation of aver- 
age raw signal level is seen, much of it following a diurnal 
pattern. The diurnal variation may be explained by the fact 
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Fig. 2 — A cross section of a typical star gives an indication of the sharpness 
and focus of the images. No change in focus or image sharpness was de- 
tected during the 10-day run. 

that only the CCD analog-clock drivers and the preamp were 
in the temperature-controlled room with the dewar and CCD. 
The ADC and other electronics were subject to daily tem- 
perature variations of 2 to 5 °F. 

4. DATA REDUCTION 

A square 28X28 pixel aperture was used to extract the 
star data. A larger aperture was not used due to crowding of 
the artificial stars. No background subtraction was performed 
because most of the background light was actually scattered 
light from the stars themselves. (Background amounted to 
about 0.5% of the total signal.) The total 1 a fractional error 
(photon noise plus detector noise) was 5.0X10 3 4 5 when av- 
eraged over the full data set after the data were normalized 
using the average signal level of all the stars. Part of this 
error is due to photon noise, which is at the level of 
1.18X 10~ 5 . We suspected that additional error would be due 
to small motions of the artificial stars across the nonuni- 
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Fig. 3— The absolute signal level seen during the 10-day run varied by about 
0.5%. Causes of variation include mechanical motions and temperature 
variations due to regular refilling of the liquid-nitrogen dewar, as well as 
some diurnal temperature variations that affected power supplies, amplifiers, 
and the analog-to-digital converter. 
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Fig. 4 — After dividing the signal from “star” 1 1 by the average signal from 
all 21 tested “stars,” the signal variation is reduced to roughly five parts in 
10 5 . Note the improved stability during the second half of the run. This is 
believed mainly due to a drop in excess dark current during the first few 
days. 

formly sensitive face of the CCD (Jorden et al. 1994), which 
would lead to reproducible errors that could be measured and 
corrected. To assess this possibility, the star images were 
individually centroided in each recorded image using mar- 
ginal analysis (Elliot et al. 1989). 

It is apparent that changes in the normalized brightness of 
a star are correlated to some extent with motions of its cen- 
troid. The periodic spike features seen in Fig. 4 correspond to 
the times when the dewar was refilled. Changes in the weight 
of the dewar, as well as local temperature changes due to the 
LN fill process, cause mechanical motion of the CCD and of 
the optical assembly. The resulting image motions [see Fig. 
5(a) for the motion of star 11] provide a well-sampled func- 
tion of brightness versus centroid position. As expected, the 
polysilicon electrode structure on the front side of the CCD 
caused a well-defined correlation (or anticorrelation) be- 
tween brightness and image position to be seen for most of 
the stars. This relationship is plotted in Fig. 5(b) for star 11, 
and it is seen that for this star, the X position is correlated 
with brightness, while the Y position is anticorrelated. In 
general, X and Y are either correlated or anticorrelated, de- 
pending on the image position with respect to the CCD struc- 
ture. For the sample of 21 stars that were analyzed, 
correlation/anticorrelation occurred roughly an equal number 
of times for X and T. Of the four stars that showed little or 
no XI Y correlation, one had the lowest uncorrected fractional 
error (1.7X10 5 ), and the remaining three showed a weak 
dependence on raw brightness, probably due to nonlinear 
system response. Figure 6 shows star 43 ’s behavior. 

Because of the strong linear correlation/anticorrelation 
observed between normalized brightness and image position, 
as well as the observed weak dependence of normalized 
brightness on raw brightness, the normalized brightness data 
for each star were fit independently with a multiple- 
parameter linear least-squares fit: 

Bfo=C y +C 2 8X+Ci8Y+C A SB. 

Here, Z? fil is the fitted value of the normalized brightness 



Fig. 5 — (a) The centroids of the “stars” moved slightly on the CCD. The 
motions were somewhat repeatable; the main motion was closely synchro- 
nized with the periodic refilling of the dewar. (b) There is a clear correlation 
between the position of the “star” (here star 11) on the CCD and the nor- 
malized brightness. This is almost certainly due to the nonuniform sensitiv- 
ity across each pixel, caused by light attenuation of the polysilicon elec- 
trodes that cover the CCD. 
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Fig. 6 — There is a small correlation between the actual raw signal level and 
the normalized brightness. This is believed due to small nonlinearities in the 
electronic system. 
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Fig. 7 — By subtracting measured corrections due to changes in signal level 
and to motion of the image on the CCD, a considerably improved residual 
error is found, as shown by the lower curve. The upper curve shows uncor- 
rected normalized data, and the middle curve shows the least-squares linear 
fit for this star (star number 7). 

of the image in question, the C,s are fitted constants and SX, 
SY, and SB are the differences between the instantaneous 
values of the X position, Y position, and Brightness for a 
given star image minus the corresponding average value 
taken over the whole data set for that star. 

The resulting fit for star 7 is plotted as the middle line in 
Fig. 7 with the normalized data plotted as the upper line for 
comparison. The difference between the data and the fit is 
plotted as the lower line in the figure. The fitting process 
provides about a factor of 3 improvement over the raw data, 
yielding a 1 a fractional error of 1.8X10 -5 and does a good 
job of flattening the data and reducing the spikes caused by 
motion that occurred when the dewar was refilled, although 
the standard deviation is still above the shot-noise limit of 
1.18X10” 5 . Note however that the fractional error improved 
with time. The error in the corrected data for the first 120 
hours of the test was 1.96X 1(T 5 and dropped to 1 .61 X 10~ 5 
for the last 100 hours. The implication is that as the system 
stabilizes with time, we approach the shot-noise limit. As- 
suming that the data from the second half of the experiment 
are representative of a stable system, then by subtracting in 
quadrature the photon-noise fractional error of 1.18X10 
(as measured for this experiment) from the total measured 
fractional error of 1.6X10 5 , we find the fractional error 
contribution from the CCD plus any uncorrected experimen- 
tal errors to be 1.1 X10~ 5 . The result of fitting all stars over 
the whole duration of the experiment is shown in Fig. 8. 

In order to determine the expected residual detector- 
related noise for the proposed long-term planetary detection 
observing program mentioned earlier, we tried using only the 
data from the last half of the test run, quadratically sub- 
tracted the photon-noise contribution, and averaged the data 
to an integration time of 5 hr. The result of this analysis is 
shown in Fig. 9. The 1 cr fractional error for all the stars in 
this case is 0.99X 10 \ 

The larger errors during the first part of the experiment 
are not surprising, since the CCD was severely overexposed 
to light while cold, during the initial setup. This produced 


Fig. 8 — The standard deviation of the corrected signal level for the 10-day 
run is approximately 1.8X10 \ Larger deviation of a few “stars may be 
due to local imperfections of the CCD. This error includes 1.18X 10 3 due 
to photon shot noise, and is at an integration time of 35 min. 

excess dark current whose decay could readily be detected 
during the first few days. The overflow and spillage of liquid 
nitrogen after the first day of operation may have also con- 
tributed uncorrected errors. 

5. DISCUSSION 

This experiment was designed to demonstrate the poten- 
tial of CCDs for a specific photometric project where very 
small and relatively short-term variations in brightness are to 
be measured for a large number of stars, and the average 
brightness of all the stars can be used to remove most 
instrument-induced variations. Absolute photometry is not an 
issue, and the background level will be very small and quite 
constant. Since each star will illuminate its own set of pixels, 
nonuniform color response should not interfere with the 
measurements. Changes lasting a few hours are of major 
interest for the planned observations. 

There may be some concern about the possible effects of 
cosmic-ray induced signals, as well as gradually increasing 
radiation damage to the CCD. However, the rate of cosmic 
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Fig. 9— The standard deviation can be further improved by using only the 
more stable data from the last half of the run, and by averaging over a 5-hr 
interval, which is more representative of the planned experiment. 
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rays in space should be about 1 s -1 cm -2 . If each star covers 
15 pixels, each 15 microns square, an area of about 3X 10~ 5 
cm 2 , there would only be two or three hits per star per day, 
an effect very much smaller than shot-noise. Radiation dam- 
age may change the charge-transfer efficiency, and may also 
lead to increased dark signal, but these changes should easily 
be detected and corrected. Of course, permanent changes in 
local sensitivity due to radiation damage may be indistin- 
guishable from permanent changes in the brightness of a star. 

We recognize that there are many other limitations to pho- 
tometric accuracy with CCDs, particularly for ground-based 
observations (see Young et al. 1991). While the data from 
this experiment demonstrate good intrinsic long-term stabil- 
ity for a front-illuminated CCD used for relative photometry 
at one color with a number of stars on the CCD, the experi- 
ment does not address numerous other problems seen in as- 
tronomical photometry. 

One noteworthy problem is that the spatial response of 
most CCDs is wavelength dependent. This is well known for 
back-side-illuminated CCDs, where the UV and blue re- 
sponse can vary dramatically over the surface of a CCD, and 
also may change over a short time for CCDs whose UV 
response is obtained by UV flooding. At a lesser level spatial 
color dependence has also been found to exist for front- 
illuminated CCDs, although no time dependence is expected. 
Differences of 1% or so between the relative response across 
a front-illuminated CCD at 650 and 700 nm wavelength light 


have been seen in our laboratory. For back-illuminated 
CCDs, the problem is usually much worse and, in addition, 
interference fringes with amplitudes of 10% or more are of- 
ten detectable at wavelengths about 800 nm. 

These effects mean that changes in the color of sky back- 
ground, for example, may make background subtraction less 
accurate than it should be. They also may require careful 
calibration of each detector at each wavelength of interest. 
However, for time-dependent relative photometry in space, 
the CCD appears capable of a relative precision of 1 part in 
10 5 and is as close to an ideal detector as we are likely to see. 

L.B.R. wishes to acknowledge support through Sverdrup 
Technology, NASA Grant No. NAS2- 13769. 
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